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Abstract — We consider a time-slotted two-hop wireless system 
in wliicli the sources transmit to the relays in the even time 
slots (first hop) and the relays forward the packets to the 
destinations in the odd time slots (second hop). Each source may 
connect to multiple relays in the first hop. In the presence of 
interference and without tight coordination of the relays, it is not 
clear which relays should transmit the packet. We propose four 
decentralized methods of relay selection, some based on location 
information and others based on the received signal strength 
(RSS). We provide a complete analytical characterization of these 
methods using tools from stochastic geometry. We use simulation 
results to compare these methods in terms of end-to-end success 
probability. 

I. Introduction 

We consider a two-hop wireless ad hoc network in which 
the sources are distributed randomly on the plane and each 
source has a destination at a distance i? in a random direction. 
In addition there exists a set of relays (different from the 
destinations) which assist the sources. In the first hop each 
source transmits and will be decoded by any node listening 
(relays and destinations) with signal-to-interference ratio (SIR) 
greater than a fixed threshold T. In the second hop, some of 
the relays which were able to decode information in the first 
hop, transmit. So in the first hop each packet may be received 
by many relays, hence multiple copies of the same packet may 
exist at different relays. In networks with low or no mobility 
multiple copies may be avoided by setting the routing tables 
apriori and relays rejecting packets if the source-destination 
pair is not in its routing table. Without routing tables, the 
relay nodes having the same packet may have to coordinate 
with each other and then choose one among them to transmit 
the packet. In a mobile wireless network such a coordination 
incurs significant overhead and also restrict the number of 
relays per source-destination pair to one. More importantly 
such a restriction may in effect reduce the probability of 
success. It is not clear how to choose a subset of these 
intermediate relays in a distributed fashion so as to reduce the 
interference and increase the probability of packet delivery. In 
this paper, we analyze the success probability of such a two- 
hop scheme taking the interference and the spatial statistics of 
the transmitting nodes into account. 

Related Work: In [3], relay selection is based on the SINR. 
All the relay-destination channel states are assumed to be 
known at the destination, and the destination chooses the 



relay. In [5] the relay selection is based on the channel state 
information that is fed back to the source. In [2] the relays 
estimate the channel using channel reciprocity theorem and 
use timers to select the best relay. In [8] a relay selection 
method (GeRaF) based on the distance from the destination 
is considered. TDMA based contention is used to resolve the 
relays transmitting to the same destination. In all these method 
some form of channel contention and feedback is used to select 
the relays. In the methods we analyze relays are selected in a 
completely distributed fashion without any channel feedback. 
This makes these relay selection schemes suitable in scenarios 
with moderate to high mobility. 

The main contribution of this paper is the complete analyti- 
cal description of a two-hop wireless network with interference 
when the nodes and relays are distributed as a Poisson point 
process on the plane. The techniques developed in this paper 
can be used to analyze other position based relay selection 
methods. 

II. System Model 

We assume that the sources form a Poisson point process 
(PPP) ips of intensity As on the plane. The relays are also 
assumed to form a PPP (pr of intensity on the plane. Each 
source x ^ (ps has a destination denoted by r{x) (not a part 
of (ps or (pr) at a distance R in some random direction. We 
assume that the fading between any two nodes is Rayleigh 
distributed so that the powers are exponential with unit mean. 
A transmitter located at x can communicate with a receiver at 
y if SIR(a;, y \ (p) > T. The SIR is defined as 

cTu/ I j,\ Ky9{x-y) 
SIR(a;, y\<p) = ^ "—r r 

where (p is the transmitting set, g{x) is the path loss function, 
and hxz is the power fading coefficient between nodes x and 
z. We assume T > 1, i.e., a narrow band system which implies 
at most one transmitter can connect to any receiver The path 
loss function g{x) is assumed to depend only on ||a;||, to 
monotonically decrease with and linia^oo — 

to guarantee finite mean interference We restrict the number 
of hops between any source-destination pair to be two. So a 
source can reach its destination either in a single hop or can 
use the relays to reach the destination. We can assume that 
the sources transmit in the even time slots and a subset of the 
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relays in the odd time slots. 

Notation: 

• We define 

l{x^v\<P) = l{SlR{x,y\cp)>T). 

l{x ^ y \ (j)) represents a random variable that is equal 
to one if a transmitter at x is able to cormect to a receiver 
y when the transmitting set is (j). 

• We define for x £ 

(i>r{x) = {y e(pr ■■ l{x ^y\ (j>s)} ■ 

(t>r{x) denotes the cluster of relays to which the source 

X is able to connect in the first hop. 
Metric: We analyze the success probabiUty for the direct 
connection between the source-destination and the two-hop 
connection between them separately. Let P\ denote the prob- 
abiUty that a source can connect to its destination directly in 
the first hop. More precisely we define 

Pi = lim -^E V l{x ^ r{x) \ cj),). (1) 

The relays which were able to connect to some source in the 
first hop are the potential transmitters in the second hop. In the 
relay selection methods studied in the next section, a subset of 
these potential transmitters A^^ C (jjrix) are selected for each 
X € 4>s transmit in the next hop. Let the probability that a 
relay can connect to its intended destination (determined by 
the source to which it cormects in the first hop) in the second 
hop be P2, defined as 

P2= lim > ; l-f{x), (2) 



lim ^^E y 1 

x^4>sC\B(o,a) 



where f{x) = IlyeAr, " ^ ''(2^) I ^)) ^nd ip = 
Ua;g0^iVj;. In the above equation 1 — f{x) is equal to one 
if and only if at least one node belonging to (j)r{x) is able 
to connect to the destination r{x). Here we are assuming no 
cooperative communication between nodes which have the 
same information, so relays belonging to the same cluster 
(l>r{x) also interfere with each other in the second hop. Then 
the success probability is given by 

= 1-(1-Pl)(l-P2). 

In the above equation we assumed that the success probability 
of the direct connection is independent of the two hop success 
probability. We also used the spatial ergodic property of PPP 
in defining Pi and P2. 

III. Location-Unaware Relay Selection 

In all the four methods described below, a relay has to make 
a decision whether to transmit in the second hop if it is able 
to connect to some source in the first hop. More precisely: 

1) The destination r{x) can directly decode a packet from 
X e (ps in the first hop if the SIR(a;, r{x) | T. 

2) A subset of the relays C (^^(a;) are chosen in a 
distributed manner to transmit in the second hop. 



3) In the second hop, the destination r{x) can decode the 
packet from a relay y € (pr ii the SIR(r(a;), y \ tp) > T. 
The success probabihties of all these methods is analyzed in 
the Appendix. 

A. Method 1: All relays transmit 

This is the most basic scheme where all relays which receive 
in the first time slot transmit in the second hop: 

Nx = (prix). 

As we shall see, this method has the worst performance 
because of the high interference present in the second hop 
(specially when the relay density is high). 

B. Method 2: RSS-based selection 

When a relay node is able to connect to a source node, the 
relay node has information about the RSS and could use that 
information to make a decision to transmit in the next hop. 
In Method 1, we do not utilize any information regarding the 
SIR received at the relay. In this method we utihze the RSS 
information to make the decision. We have RSS = S+I where 
S is the strength of the desired signal and I is the interference 
observed. Since the relay was able to decode the source, we 
have f > T and thus 

RSS 



I < 



S 

1 + T -f' 



Hence a small value of implies low interference (and 
hence may see few interferer's in the second hop). This might 
also mean smaller S which implies that the relay is far from 



the source. A large value of implies a large S and hence 
would indicate a relay close to the source. This may potentially 
also imply more interference at the relay. So in the second hop 
we give higher priority to nodes that observe smaller value of 
^p|t. a relay y e 4>r{x) transmits the packet in the second 



time slot with probability 



exp 



) l£i 

1 + T 



(3) 



where RSS(j/) is the RSS that the relay y observes and 6 
represents a parameter to be chosen so that 



SRSSjy) 
1+T 



IS not 



too large. We have chosen an exponential penalty just for 
convenience and the effect of the penalty function should be 
investigated further. So we have 

.RSS{y) 



Nx = <.y & (prix) ■■ Uy < exp 



1+T 



(4) 



where Uy is a set of i.i.d uniform random variables in [0, 1]. 

We observe that Methods 1 and 2 are completely decen- 
tralized and require no information about the location of any 
node. Hence these methods of relay selection work perfectly 
even with high mobility and incur zero overhead. 

IV. Location-aware selection 

In the location-aware based methods we assume that each 
node has knowledge about its own location and each source 
knows the spatial location of its destination. Also each packet 
has information about the location of the source from which 
it originated and the location of its destination in its header. 
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A. Method 3; Sectorized relay selection 

After a relay receives a packet in the first hop, it calculates 
the angle between the relay-source and the source destination 
and makes a decision based on this information. More pre- 
cisely a relay y £ (prix) transmits the packet in the second 
time slot if Zyxr{x) < 9. So for this method, we have 

In this method we are reducing interference by choosing relays 
in a sector If the angle 9 is properly chosen the sector may 
contain only one relay, i.e., \Nx\ = 1 which would eliminate 
the intra-cluster interference. 



B. Method 4 .• Distance -based selection 

In this method, we select the relays in (f>r (x) depending on 
their distance from the destination r{x). A relay y G (f>r{x) 
transmits the packet in the second time slot with probability 

2\\y-r{x) 



X =1,g(x)=|xr'',e=7t/6 



exp 



R 



(5) 



We have normalized the distance by R/2 so that nodes closer 
to the source are not over penalized, e is a normalizing 
parameter which we will choose later. So we have 

N^ = [ye ^rix) : Uy < exp (-e2||y - r{x)\\R-^)} 

where Uy is a set of i.i.d uniform random variables in [0, 1]. 
In this method, we observe that relays close to the destination 
have a higher chance of transmitting than those closer to the 
source. One could replace ||y — r(a:)|| in (|5]l with \\y — x\\ 
and exp with 1 — exp. Then the source need not know its 
destination location but one would loose the directionality of 
relay selection. 

In the location-aware methods proposed in this section, 
each source needs to know where its destination is located 
and maintaining this information would require a significant 
overhead in a mobile network. 
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Fig. 1. Comparison of Method 3 versus tlie center-sclieme. Tlie success 
probability versus the relay intensity Ar is plotted for various values of the 
source destination distance R. 
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Fig. 2. Comparison of the simulation and theoretical results for Method 3. 



V. Simulation Results 

In this section we compare the four methods described in the 
above section by simulations. For the purpose of simulation 
we consider a square [—30, 30]^ in which the nodes are located 
and use Monte-Carlo method to evaluate the results. We use 
T = 3 for all the simulations. 

In Figure [T] we compare Method 3 with a scheme in which 
the relays are apriori chosen. In this scheme each source- 
destination pair has one relay assisting in communication. 
The relay is centered halfway between the source and the 
destination. Intuitively such a scheme is the best single-relay 
scheme in terms of end-to-end success probability. We observe 
that using uncoordinated relays yields a better performance 
than selecting a relay apriori. This is because of the selection 
diversity that occurs due to fading and the node locations. 
Also observe that there is an optimal value of relay intensity 
that achieves the maximum value of P,. In practice this can 
be achieved by starting with large number of relays and 
using ALOHA-like thinning in the second hop. We observe 



that the intensity of relays required to outperform the center- 
method decreases with increasing source destination distance 
R. From Figure |2] we observe that the probability of success 
as computed from the theory matches the simulations closely. 
In Figure [3] the success probability is plotted for the four 
methods described in previous sections. We first observe that 
all the methods have an optimal value of the parameters 8, 
9 and e that achieve the maximum success probability. We 
also observe that Method 1 has the worst performance. This 
is because of the high interference caused by all the relays 
(that have received a packet in the first hop) transmitting in 
the second hop. We observe that Method 3 i.e., the sectorized 
selection method performs the best. Also observe that the RSS 
based selection has twice the Pg as compared to Method 1 in 
which every relay transmits. We also observe that the both the 
location-aware schemes outperform Method 1 and Method 2 
for the particular values of T = 3, — 1, A,. = 1.5. 



4 



R=1 .4, g(x)=|xr , T=3,). =1 , X =1 .5 
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Fig. 3. Probability of success for the four methods described in the previous 
section. The x-axis represents 5 for Method 2, A8/n for Method 3 and e/6 
for Method 4. 



VI. Conclusion 

In this paper, we have have analyzed the success probability 
in a two hop system taking interference and spatial distribution 
of the nodes into account. We have provided an analytical 
solution for all the four methods using some approximations. 
This method of analysis can be easily extended to any position 
based relay selection. We have shown that uncoordinated se- 
lection of relays increases the success probability as compared 
to selecting a relay for each source-destination pair apriori. 

Appendix 

The probability that a destination located at z can decode 
the packet transmitted by a source £, when the interference is 
caused by (j)s is 



- p(siR(^,z|0.)>r) 

= ^^hz > — ^^z9{x - z)) 



(6) 
(7) 



(a) 



(b) 



E cxp 



9{S. - -A 



T 



5(C - z) 



hxzgix - z) 



where 



exp{-X, / I3{i-z,y)dy) 

B(x,y) = 



(a) follows from the exponential distribution of h^^, and (b) 
follows the probability generating functional (PGFL) of the 
PPP [1], [7]. Also observe thatp^(z) depends only on ||^ — 
Direct transmission: So from ^ we have 

Pi = lim y 1(.T ^ r{x) \ 

xe<l)snB(o,a) 

= lim ^— / px{r{x))dx = Po{R), 

Q-oo TTfl^ JB{o,a) 



where (a) follows from the Campbell-Mecke theorem [7] and 
(6) follows from he fact that p^{z) depends only on ||^ ^ z\\. 

First hop: A point process is completely characterized by 
its PGFL and so we will evaluate the PGFL of the relays 
which can connect to source ^ G 0s, i-C-, the cluster <j)r{£.)- 
Let < ^(a;) < 1. The PGFL of 4)riO is given by 

GM^)) = E II 1 - (1 - vix))li^ ^ X \ (f>s) (8) 



(b) 



cxp —A,. / (1 — v{x))l{S, — > X I (f>s)dx 



> exp ~Xr / (1 — v{x))p^{x)dx 



(9) 



where (a) follows since (f>r is a PPP, and (6) follows from 
Jensen's inequality. From the PGFL we observe that the point 
process consisting of relays which connect to the origin is not 
a PPP. (6) would have been an equality if 1(^ ^ a; | (ps) are 
independent for different x and the resulting process would 
be a PPP. But for the sake of analysis, we make the following 
assumptions and justify them by simulations. 

1) We assume that the spatial distribution of 4>r{£,) is an 
inhomogeneous PPP with intensity XrP^{£,)- Since T > 

2) We also assume <f>r{£,i) is independent of (i>r{£,2) for all 

We will show the results obtained by this assumption are close 
to the actual by simulation. From Figure |2] we observe that 
the simulation results (for Method 3) are very close to that 
predicted by theory making the above assumptions. This is 
intuitive since many terms in (O are independent and thus the 
bound in (|9|l is very tight. 

A subset of relays C 0^(6 for each ^ G 0s transmit 
in the second hop depending on the relay selection method. 
This is basically a thinning of the point process (t>r{C)- We will 
now derive the intensity of the point process for different 
methods. We will denote the spatial intensity of by Aj(z) 
and we have ¥.[\N^ n A\] = /\^{z)dz for any AdR^ . 
Method 1: Since iV^ — 4>r{£,)' we have A^{z) = XrP^{z). 
Method 2: From (fTTT l. we have 



A^(z) = A.E (^exp (^-<5^^) 1 ^ z I 0s) 

{_X r {5g{i~z)+T)g{y) , 

V Jr2 g(5_2) + (5g(5_^)+T)g(a) «y 



A^(i + r)- 



exp I 



l + T + 5giC~z) 



where (a) follows from a procedure similar to the evaluation 
of P^z). 

Method 3: Given ^ and r{^), we have 

Ai{z)^XrPdz)lUz^r{0<e). 



Method 4: Given ^ and r(6, we have 

Ae(z) - XrPdz) exp {-e\\z - r(e)||i?-i) 

The average number of relays in a cluster that transmit in 
the second hop is J A^{z)dz. 
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Second hop: The transmitting set in second hop is given 
by V' = U^e0 -^i- Since T > 1, at most one transmitter 
belonging to iV^ can connect to r(^). So the probability that 
no node from A'^ can connect to r{S) denoted by /(^) is given 
by 



1 



1 {z ^ r(e) I ^) . 



So P2 (same as P2 without the limit) is given by 
P2 



E ^ 2^ 1 -> r{0 I V') 

x&<ptC\B{o,a) z£Nc: 



(a) 



Xt [ E ^ l(z-.r(C) I A^cUV')rfe 

JB{o,a) ^g^^ 

Xt f [e [A^{z)1 [z r(0 I iVf U V')] dzd^ 

JB(o,a) J 



where (a) and (b) follow from the Campbell-Mecke theorem 
and Slivnyak's theorem. We included Aj(x) in the expectation 
operator because in Methods 3 and 4, A^(z) depends on the 
random variable r(^). We now show that the inner integral 
does not depend on ^. We also have A^(z) = Ao{£, — z) for 
Method 1 and 2. For Methods 3 and 4 we have A^(z) = 
Ao(C — z) where the equality is in distribution. Using the 
substitution z' ^ z — ^, the stationarity of -0, and the above 
property of A^(z) we have 

^2 = 7^/ / E [Ao(z)l (z ^ i?^ I V U TVo)] dzdz/ 

(10) 

where i?^ = (-Rcos(i/), i?sin(z/)). For Methods 1 and 2 by 
the isotropic nature of No and we have 



Ao(z)El {z ^ R\'^l)y^No)dz 



(11) 



Due to space constraints we will only describe how to derive 
P2 for Method 3. Method 4 can be analyzed in a similar 
fashion. From ( fTOl i and the definition of Ao(z) for Method 
3, we have P2 = 

— / / Po{z)El{AzoR^ < 61)1 (z i?^ I V U iVo) d^d;/ 
27r Jo J 

Since -0 is isotropic we have 

P2 = ^ j Poiz)El -> i? I V U No^ dz. (12) 

Here No is equal to A^o n S{o, R, 6) where S{o, R, 9) denotes 
a sector of angle 9 on either side of the line joining the origin 
and (i?, 0). With a shght abuse of notation we will denote No 
also by No and No depends on the relay selection method. 
We now evaluate E [1 (z — » i? | ?/;')] where 'ijj' — -tj^ij No- 



= ¥{h,B,g{z- R)>TIiiP\R)) 
Ee.J-Il^) (13) 



E [1 (z ^ i? I V')] 



.g(z - i?) ; 

where (a) follows from the exponential distribution of h^R- 
Also ( fT3] l is the Laplace transform of the interference evaluated 
at T/g{z — R). By our assumptions ip' is a Poisson cluster 



process [4], [6] with an additional cluster at the origin. The 
Laplace transform of the interference in this case is given by 
Ecxp (-sJ2y^^' hyRgiy - R)) which is equal to 



E 



Yl exp{~shyRg{y - R)) 



1 + sg{y - R) 



where G^(.) is the PGFL of the process ip' and (a) follows 
by Laplace transform of the fading. So we have 



;[i(z ^ i? I V') 



(a) _ 

= E 



E 



n n 



, Tg(y-R) 
€60. ydN^ ^ ^ g(z-R) 



n exp (-/3(z-i?,0 



n 



, Tg(y-R) 
yeN„ ^ g{z-R) 



exp 



(-/3(z-i?,i?)) 



exp ( -(3{z - R,R)-Xs 1- exp (-/3(z - R, 



where (a) follows from assumption 2, (6) follows from a 
technique similar to the derivation of p^{z), (c) follows from 
the PGFl of PPP and where 



R2 



/3{z,y)Ao{y + 0dy. 



From the above equation, (fT2] i and (fTTT i we have, 



^2 = / A„(z + i?) 

iR2 



■ exp ( -^(z, R)-Xs J ^l-exp (-/3(z, C)) 



dz, 



where Ao(z) = Ao(z) for Method 1 and 2. For Method 3, 

Ao(z) = ^Xrp^iz). For Method 4 

Ao(2) = :fpd^) / exp (-e||z - R,\\R-^) dv. 
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